1 Introduction

The monarch butterfly (Danaus plexippus) has long captured both scientific and public attention as one of the most charismatic insect species. Its extraordinary, continent-scale migration across North America is unique among butterflies, spanning up to 4,000 km and involves multiple successive generations that complete different stages of the annual cycle (Brower, 1995). This migration is tightly coupled to reproductive diapause and seasonal timing, forming one of the most complex life-history strategies known in insects (Agrawal et al., 2014; 2021; 2025). Equally striking is the monarch’s specialized larval dependence on milkweeds (Asclepias spp.), which provide not only nutrition but also chemical protection. Monarch caterpillars ingest toxic cardenolides and sequester them through metamorphosis, rendering both larvae and adults unpalatable to vertebrate predators-a classical example of coevolutionary adaptation and aposematism (Agrawal et al., 2012; 2024; 2025).

For much of the twentieth century, research on monarchs was rooted in natural history, ecology, and physiology - from field studies on migratory navigation (Brower, 1995) and overwintering colonies in Mexico (Rendón-Salinas et al., 2023) to biochemical investigations of cardenolide sequestration (Agrawal et al., 2021; 2024). These ecological and behavioral foundations laid the groundwork for monarchs to emerge as a model system in evolutionary biology (Brower, 1995; Oberhauser et al., 2015). The genomic era began with the publication of the first monarch draft genome (Zhan et al., 2011), which provided a platform for identifying genes and regulatory networks underlying migration, chemical defense, wing patterning, and reproductive physiology. Since then, monarch research has become increasingly integrative, bridging genomics, neuroscience, chemical ecology, and conservation biology.

Advances in genomic and molecular tools have transformed monarch research, enabling detailed investigation of the genetic and regulatory bases of complex traits (Li et al., 2025). High-quality reference genomes, long-read sequencing, and chromosome-scale assemblies now provide near-complete maps of coding and non-coding regions, structural variants, and neo-sex chromosomes (Zhan et al., 2011; Mongue et al., 2017; Zhan et al., 2020). Transcriptomic approaches, including tissue-specific and developmental-stage RNA sequencing, reveal dynamic gene expression patterns underlying migration, diapause, and chemical defense. Functional genomics techniques such as RNA interference, TALENs, and CRISPR/Cas9 genome editing have made it possible to experimentally validate candidate genes and regulatory elements, linking genetic variation to phenotypic outcomes (Markert et al., 2016; Zhang and Reed, 2016). Complementary population genomic analyses, leveraging resequencing, Genome Wide Association Studies (GWAS; Uffelmann et al., 2021), and selection scans, are uncovering polygenic architectures for migratory behavior, adaptation to host-plant chemistry, and responses to parasites. Together, these tools provide a comprehensive framework for connecting genotype to phenotype, enabling integrative studies that span molecular mechanisms, ecological interactions, and conservation applications (Li et al., 2025).

In this review, we synthesize current knowledge of monarch ecology, genetics, and metabolomics (Figure 1) across seven major domains: (1) the development of genomic resources and functional tools, (2) the genetic basis of migration and seasonal behavior, (3) molecular evolution of chemical defense and host-plant specialization, (4) sex-chromosome and structural genome evolution, (5) eco-genomic interactions involving parasites, microbiota, and host plants, (6) population genomic signals that inform conservation, and (7) metabolomic perspectives on host-plant chemistry and monarch sequestration. We also highlight methodological frontiers - from long-read pan-genomics to CRISPR-based functional validation, single-cell genomics, and integrated metabolomics - that are poised to resolve outstanding questions about how the monarch’s remarkable life history is encoded in its genome and shaped by ecological interactions. Overall, this review aims not only to synthesize recent advances in monarch genomics and chemical ecology, but also to identify conceptual and methodological gaps that must be addressed to link molecular mechanisms to ecological function and conservation outcomes.

2 Genomic Resources and Technical Toolkit

2.1 Reference genomes, assemblies, and databases

The initial monarch genome assembly (~273 Mb, ~16 866 predicted protein-coding genes; Table 1) provided a foundational resource for trait mapping and comparative genomics in Lepidoptera (Zhan et al., 2011). Subsequent chromosome-level assemblies, improved annotations, and curated community databases have substantially expanded the utility of monarch genomics for population, functional, and evolutionary analyses (MonarchBase team, 2012; Mongue et al., 2017). These resources have enabled identification of candidate loci underlying migration, circadian regulation, detoxification, pigmentation, and host-plant interactions, as well as comparative analyses of sex chromosome evolution within Danaus.

Recent high-quality assemblies have further resolved structural variants, repetitive regions, and the neo-Z chromosome, providing insight into genomic features that may contribute to adaptation and phenotypic divergence among migratory and non-migratory populations (Table 1). Collectively, these monarch-specific genomic resources now support both hypothesis-driven functional studies and broad-scale evolutionary inference.

2.2 Functional genomics, genome editing, and multi-omic tools

Functional genomic studies in monarchs have leveraged transcriptomic, proteomic, and metabolomic data to link genetic variation with key migratory, physiological, and defensive traits (Table 1). For example, tissue-specific transcriptomic analyses of antennae revealed circadian clock gene expression patterns essential for time-compensated sun-compass navigation (Merlin et al., 2009), while expression profiling of fat body and flight muscle tissues clarified metabolic shifts associated with long-distance migration and lipid utilization (de Roode et al., 2011).

Monarch research has also been at the forefront of integrating chemical ecology with genomics. Metabolomic analyses of larvae and adults have quantified cardenolide sequestration and detoxification, revealing how host-plant chemistry interacts with monarch genotype to shape toxin resistance and performance (Agrawal et al., 2012; Petschenka et al., 2013; Dreisbach et al., 2023). Population genomic resequencing across migratory and resident populations has identified loci associated with neural, metabolic, and endocrine function, providing candidates for functional validation (Zhan et al., 2014).

Genome editing approaches, including TALENs and CRISPR/Cas9, have enabled direct tests of gene function in monarchs and related Lepidoptera, allowing researchers to move from association-based inference toward causal understanding of traits such as circadian regulation and pigmentation (Markert et al., 2016; Zhang and Reed, 2016). Together, these monarch-focused functional and multi-omic applications illustrate how established technologies can be combined to address uniquely integrative questions spanning behavior, chemistry, and ecology.

3 Genetic Architecture of Migration and Seasonal Behavior

3.1 Migratory phenotype complexity

Monarch migration is a compound phenotype encompassing sun-compass orientation, reproductive diapause, lipid storage, flight endurance, and seasonal timing (Table 2). Sun-compass orientation depends on antennal circadian clocks and neural integration of solar cues, as demonstrated by experimental disruption of antennal function or clock genes (Merlin et al., 2009; 2020). Reproductive diapause is regulated by photoperiod and hormone signaling (Fleming and Alto, 2006; Green et al., 2019; Freedman et al., 2023), while lipid storage and flight endurance are shaped by metabolic and mitochondrial pathways associated with long-distance flight capacity (de Roode et al., 2011; Zhan et al., 2014). Seasonal timing of migration and breeding integrates environmental cues such as day length and temperature across generations. Collectively, these traits arise from interacting circadian, endocrine, metabolic, neural, and developmental modules, reflecting a modular and polygenic architecture underlying the monarch’s migratory phenotype.

3.2 Antennal clocks, sun-compass orientation, and molecular candidates

Peripheral circadian clocks located in the antennae are essential for time-compensated sun-compass navigation, allowing monarchs to adjust orientation as the solar position changes throughout the day (Merlin et al., 2009; Hemstrom et al., 2025). These clocks interact with central brain regions involved in spatial and sensory integration, including the central complex and optic lobes (Merlin et al., 2009; Guerra et al., 2012). Core clock genes (e.g., period, timeless, cryptochromes) coordinate rhythmic gene expression underlying behavioral timing, while transcriptomic analyses implicate broader networks of sensory and metabolic genes contributing to energy allocation and flight performance (Zhan et al., 2011; Zhan and Reppert, 2013). Together, these findings highlight migration as an emergent property of interconnected circadian, sensory, and metabolic pathways rather than a single master regulator.

3.3 Population genomics: migratory vs. resident populations

Comparative genomic analyses of migratory and resident monarch populations reveal modest but consistent differentiation at loci associated with neural signaling, lipid metabolism, and endocrine function, consistent with adaptation to migratory versus sedentary lifestyles (Zhan et al., 2014; de Roode et al., 2013). Divergence at endocrine-related loci involved in juvenile hormone signaling may contribute to differences in reproductive diapause between populations (Hemstrom et al., 2025). Most signals are dispersed across the genome, supporting a largely polygenic architecture for migration, although a small number of loci with larger effect sizes may disproportionately influence key migratory traits (Freedman and Kronforst, 2023).

3.4 Toward causal tests

Advances in genomic mapping, selection experiments, and functional genetic approaches provide promising avenues for testing causal links between genotype and migratory phenotypes. Genome-wide association studies, artificial selection lines, and pedigreed crosses can be combined with CRISPR/Cas9-based validation to interrogate candidate genes involved in circadian regulation, sensory processing, and energy metabolism (Di Cristina et al., 2025). Successful genome editing in monarchs and related Lepidoptera demonstrates the feasibility of moving from genotype-phenotype correlations toward mechanistic understanding of migratory and sensory traits (Markert et al., 2016; Zhang and Reed, 2016).

4 Chemical Defense: Na+/K+-ATPase Evolution and Cardenolides

Monarch larvae sequester cardenolides from milkweeds, which bind and inhibit Na⁺/K⁺-ATPase (ATPα; Mongue et al., 2025). Specific amino-acid substitutions in ATPα, such as N122H and Q111L, reduce binding affinity for cardenolides and confer resistance (Petschenka et al., 2013; López-Goldar et al., 2024). Convergent evolution has been observed in other specialist herbivores (Agrawal et al., 2024), like Danaus chrysippus and Tetraopes beetles, which carry similar substitutions conferring toxin resistance. Biochemical assays have demonstrated that these substitutions maintain ion pump function while reducing cardenolide binding, illustrating a clear genotype-phenotype link. Variation in ATPα selectivity among monarch populations correlates with milkweed species in their breeding ranges, highlighting an eco-genomic interaction between host-plant chemistry and monarch defense strategies (Petschenka et al., 2013; Agrawal et al., 2012; 2024).

Metabolomic studies of milkweed species reveal extensive variation in cardenolide composition and other secondary metabolites, which can interact with monarch genotype to influence sequestration efficiency, developmental success, and predator deterrence (Malcolm and Brower, 1989; Petschenka and Agrawal, 2015; Agrawal et al., 2024). Integrating chemical and genomic data allows a detailed understanding of coevolution between monarchs and their host plants.

5 Pigmentation, Mimicry, and Developmental Genetics

Monarch wing coloration exemplifies aposematism, with a consistent orange-black pattern across populations that reinforces predator avoidance in conjunction with cardenolide-based chemical defense. Although pigmentation genetics is not known to contribute directly to migratory behavior, it plays an important indirect role in defense by enhancing the effectiveness of warning coloration and mimicry systems. Functional studies suggest that genes such as optix, cortex, and WntA, known from Heliconius butterflies, regulate pigment deposition and patterning in monarchs (Martin et al., 2012; Livraghi et al., 2024; 2025). Regulatory changes in cis-elements of these loci likely modulate spatial expression of pigment genes rather than altering protein sequences directly, a pattern observed broadly across Lepidoptera (Ben Chehida et al., 2025). Transcriptomic analyses of wing discs have identified differential expression of yellow and dopa decarboxylase during late pupal stages, linking enzymatic pathways to melanin and ommochrome deposition (Shen, 2024). Comparative studies across Danaus species may further elucidate the molecular basis of subtle color pattern variation relevant to predator learning and mimicry (De-Kayne et al., 2025).

6 Sex Chromosomes and Structural Genome Evolution

Monarchs possess a neo-sex chromosome (neo-Z) resulting from an autosome-Z fusion approximately 5-10 million years ago (Mongue et al., 2017; Mora et al., 2024). Genes on the neo-Z show accelerated evolution and increased sex-biased expression compared to autosomes. For example, doublesex and genes involved in hormone signaling are enriched on the neo-Z (Mora et al., 2024), potentially influencing sexual dimorphism in wing coloration and reproductive timing. Long-read sequencing has identified inversions and structural variants on the neo-Z that may reduce recombination and preserve co-adapted gene complexes related to migration and diapause. Transposable element (TE) accumulation on sex chromosomes also contributes to structural diversification, shaping the evolutionary trajectory of Danaus genomes (Davey et al., 2016).

7 Parasites, Microbiome, and Eco-genomic Interactions

The specialist protozoan Ophryocystis elektroscirrha (OE) reduces monarch survival, flight performance, and fecundity (Altizer and Oberhauser, 1999; Agrawal et al., 2012). Studies show that monarchs feeding on high-cardenolide milkweeds carry lower parasite loads, linking chemical defense to disease resistance (de Roode et al., 2008; Müller-Theissen et al., 2025). Metabolomic analyses of milkweed secondary compounds, including digitoxin, asclepin, and calotropin, reveal variation in toxicity and sequestration efficiency across populations. Monarch gut microbiomes, dominated by Enterococcus and Lactobacillus species, influence digestion and detoxification pathways (Sanaei et al., 2024); RNA-seq profiling suggests microbial modulation of host immunity and nutrient assimilation (Hammer et al., 2014; van der Hoeven et al., 2013). Integration of host, parasite, plant chemistry, and microbiome data provides a holistic view of eco-genomic dynamics.

8 Population Genomics, Demography, and Conservation Relevance

Large-scale population genomic analyses indicate that eastern North American monarchs exhibit modest population structure, with subtle differentiation between migratory and non-migratory populations (Zhan et al., 2014; Freedman and Kronforst, 2023). Loci associated with lipid metabolism, circadian regulation, and neural function show signals consistent with local adaptation linked to migratory behavior. Demographic modeling further suggests that overwintering populations in Mexico have declined by more than 80% over the past two decades, underscoring ongoing conservation concern (Semmens et al., 2016; Thogmartin et al., 2017). Genomic diversity metrics are increasingly used to identify conservation-relevant units and track the retention of adaptive variation related to migration and reproduction. Integrating these genomic insights into management strategies-such as habitat restoration and protection of migratory corridors-remains an important but unresolved challenge (U.S. Fish and Wildlife Service, 2020; Erickson et al., 2023).

9 Outstanding Gaps and Priority Directions

Despite remarkable advances in monarch genomics and functional biology, several critical gaps remain that limit our mechanistic understanding of migration, chemical defense, and adaptation to environmental stressors.

9.1 Causal mapping of migratory behavior

Migration in monarchs is increasingly recognized as a polygenic and modular trait shaped by circadian regulation, neural circuitry, endocrine signaling, and metabolic pathways (Table 2). Although population genomic and quantitative genetic studies have identified candidate loci associated with migratory components, direct tests of causality remain limited (Freedman and Kronforst, 2023). Integrating high-resolution mapping with functional validation-such as allele-specific perturbations targeting genes implicated in orientation, diapause, lipid storage, and flight endurance-will be essential for resolving how combinations of alleles produce coordinated migratory phenotypes (Markert et al., 2016).

9.2 Regulatory and cell-type-specific neurogenomics

Migratory orientation and seasonal behavior are mediated by complex neural circuits, particularly in the antennae and central brain. Single-cell RNA sequencing (scRNA-seq), spatial transcriptomics, and ATAC-seq of neural tissues under different photoperiods and temperatures can reveal cell-type-specific regulatory programs. Integration with functional assays will clarify how environmental cues are encoded at the molecular level to influence migration timing and sun-compass navigation (Merlin et al., 2009; Guerra et al., 2012; Agrawal et al,. 2024).

9.3 Pan-genome and structural variation analyses

Recent discoveries of chromosomal rearrangements, including a neo-Z chromosome, highlight the potential importance of structural variation in monarch adaptation (Höök et al., 2024). However, most genomic inferences still rely on a single reference genome. Pan-genome approaches incorporating migratory and non-migratory populations will enable systematic assessment of structural polymorphisms, copy-number variation, and transposable elements that may influence migration, diapause, or chemical defense (Davey et al., 2016; Mongue et al., 2017). A key unresolved question is how such variants contribute to adaptation through regulatory versus coding effects.

9.4 Host-microbe-parasite-plant interactions

Monarch fitness and migration are tightly linked to interactions with milkweed chemistry, gut microbiota, and the protozoan parasite Ophryocystis elektroscirrha (OE; Müller-Theissen et al., 2025). Controlled factorial experiments combining milkweed species with differing cardenolide profiles, manipulations of the microbiome, and OE exposure can dissect the genetic and metabolic bases of host tolerance, detoxification, and parasite resistance (de Roode et al., 2008; Satterfield et al., 2015; Müller-Theissen et al., 2025).

9.5 Metabolomics and chemical ecology of host plants and monarchs

A particularly understudied frontier is integrating metabolomics to link plant chemistry, monarch metabolism, and stress responses. Milkweeds produce a diverse array of cardenolides, alkaloids, and other secondary metabolites (Hoogshagen et al., 2024) that vary among species, populations, and environmental conditions. Monarch larvae ingest and sequester these compounds (Agrawal et al., 2012; Betz et al., 2025), but the dynamics of sequestration, biotransformation, and excretion remain poorly characterized. Mass spectrometry-based metabolomics, including LC-MS/MS and untargeted metabolite profiling, can quantify both plant and insect chemical landscapes, revealing:

How monarchs metabolically adjust to high vs. low cardenolide diets.

The impact of environmental stressors (temperature, drought, pesticide exposure) on metabolite accumulation and detoxification pathways.

Relationships between metabolite profiles and parasite resistance, flight performance, or survival during overwintering.

Combining metabolomics with transcriptomics and proteomics can map the regulatory networks linking detoxification enzymes (e.g., cytochrome P450s, ABC transporters) to physiological and behavioral traits, providing a mechanistic understanding of how monarchs cope with chemically complex host plants.

9.6 Conservation integration

Translating molecular, genomic, and metabolomic insights into actionable conservation strategies remains a critical priority for monarch research (Box 1). Although habitat loss and climate change are widely recognized drivers of monarch declines, increasing evidence suggests that the quality, chemical composition, and genetic compatibility of restored habitats may be as important as habitat quantity (Thogmartin et al., 2017; Erickson et al., 2023). Milkweed species and genotypes vary substantially in cardenolide concentration, chemical composition, and inducibility under environmental stress, with downstream effects on monarch detoxification, growth, and resistance to the protozoan parasite Ophryocystis elektroscirrha (OE; Agrawal et al., 2012; de Roode et al., 2008; Petschenka and Agrawal, 2015). Metabolomic profiling of candidate milkweeds used in restoration programs could therefore identify plant chemotypes that balance larval performance, chemical defense, and disease suppression, while avoiding unintended consequences such as selecting for excessively toxic plants that impair development or migratory performance (Malcolm and Brower, 1989; Agrawal et al., 2024; 2025).

At the population level, coupling genomic and metabolomic data across landscapes offers a complementary framework for conservation planning. Population genomic analyses reveal consistent signals of selection on loci associated with migration, lipid metabolism, circadian rhythms, and endocrine regulation (Zhan et al., 2014; Freedman and Kronforst, 2023), suggesting that adaptive capacity may erode before demographic declines become evident. Integrating these genomic markers with metabolomic signatures-such as cardenolide sequestration profiles, lipid reserves, and stress-response metabolites-could help identify populations that retain key functional traits required for long-distance migration and overwintering success, while providing early warning indicators of declining resilience (Semmens et al., 2016; Thogmartin et al., 2017). Together, advances in causal genomics, pan-genomics, regulatory neurogenomics, and metabolomics position monarchs as a powerful model for linking genotype, chemistry, behavior, and fitness across ecological scales, informing conservation strategies that preserve both population size and the functional diversity required to sustain migration under ongoing environmental change (U.S. Fish and Wildlife Service, 2020; Erickson et al., 2023).

10 Recommended Approaches and Experimental Roadmap

To advance understanding of monarch migration, chemical defense, and adaptation, research should integrate high-resolution genomic, functional, and ecological approaches. Large, replicated mapping cohorts, including pedigreed crosses and population-resequencing panels, can identify loci underlying complex traits such as migratory orientation, diapause, and lipid storage. Functional validation via CRISPR/Cas9 or TALEN-mediated knockouts and allele swaps will allow causal testing of candidate genes. Complementary single-cell and spatial transcriptomics of antennae and brain tissues can reveal cell-type-specific regulatory programs controlling navigation, circadian timing, and sensory processing.

At the genomic level, long-read sequencing combined with Hi-C scaffolding enables chromosome-scale assemblies and pan-genome analyses, facilitating discovery of structural variants (Livraghi et al., 2024), neo-sex chromosome polymorphisms (Mongue et al., 2017), and transposable element dynamics (De-Kayne et al., 2025). SV-aware association testing can link structural changes to adaptive traits. Integrating these genomic data with metabolomic and proteomic profiling will illuminate the biochemical pathways through which monarchs process milkweed toxins (Agrawal et al., 2012; 2024; 2025) and respond to environmental stressors (Dalla et al., 2014), including parasite exposure (Altizer et al., 2015) and habitat change (Green and Kronforst, 2019).

Finally, experimental designs should bridge molecular and ecological scales. Field and laboratory studies combining genotype, microbiome composition, host-plant chemistry (Dale and Stumpe, 2014), and parasite load can quantify the ecological relevance of genetic and metabolic variation. This integrative roadmap positions future research to resolve mechanistic links from genotype to phenotype to fitness, while providing actionable insights for conservation management, such as selecting milkweed species and populations that optimize monarch survival, detoxification capacity, and migratory performance.

11 Conclusions

Research on the monarch butterfly has progressed from classical natural history and ecological observations to an integrative genomic, functional, and eco-physiological understanding of this iconic species. Chromosome-scale genome assemblies, population resequencing, and functional genomic tools have enabled identification of genes and pathways associated with migratory behavior, circadian rhythms, chemical defense, and sex-chromosome evolution (Satterfield et al., 2015). Studies of host-plant interactions and metabolomics have revealed how cardenolide sequestration and chemical stress responses mediate survival, predator avoidance, and parasite resistance, highlighting the complex interplay between genotype, phenotype, and environment (Agrawal et al., 2012; 2024; 2025). Simultaneously, ecological and population genomic analyses underscore how microbiomes, parasite pressures, and environmental changes shape adaptive variation and influence conservation priorities (de Roode et al., 2008; Dale et al., 2014; Sanaei et al., 2024).

Despite these advances, significant gaps remain in fully elucidating the proximate and ultimate mechanisms underlying monarch adaptation. Key opportunities include causal mapping of migration-related alleles, single-cell and spatial genomics to resolve cell-type specific regulatory networks, pan-genome analyses of structural variants and neo-sex chromosome evolution. Additional metabolomic profiling to link host-plant chemistry with physiological stress responses. Integrative experimental designs combining genotype, metabolome, microbiome, milkweed chemistry, and parasite exposure (Dreisbach et al., 2023; Agrawal et al., 2025) offers a pathway to connect molecular mechanisms with ecological function. By leveraging these multi-dimensional approaches, future research can provide mechanistic insights into migration, chemical defense, and adaptation while directly informing conservation strategies to enhance monarch resilience in the face of habitat loss, climate change, and shifting ecological pressures.

Acknowledgments

I thank colleagues and mentors at Kansas State University and Central College for insightful discussions that helped shape the ideas synthesized in this review. I am also grateful to the anonymous reviewer for their thoughtful and constructive feedback, which improved the clarity and scope of the manuscript. I thank the researchers in the monarch biology, genomics, and chemical ecology communities whose foundational and ongoing work made this synthesis possible. No specific funding sources directly supported the preparation of this manuscript.

Conflict of Interest Disclosure

The author declares no conflicts of interest.

References

Agrawal A.A., Petschenka G., Bingham R.A., Weber M.G., and Rasmann S., 2012, Toxic cardenolides: chemical ecology and coevolution of specialized plant-herbivore interactions, Annual Review of Ecology, Evolution, and Systematics, 43: 59-80.

Agrawal A.A., Böröczky K., Haribal M., Hastings A.P., White R.A., Jiang R.W., and Duplais C., 2021, Cardenolides, toxicity, and the costs of sequestration in the coevolutionary interaction between monarchs and milkweeds, Proceedings of the National Academy of Sciences, 118(16): e2024463118.

https://doi.org/10.1073/pnas.2024463118

Agrawal A.A., Hastings A.P., and Duplais C., 2024, Testing the selective sequestration hypothesis: Monarch butterflies preferentially sequester plant defences that are less toxic to themselves while maintaining potency to others, Ecology Letters, 27(1): e14340.

https://doi.org/10.1111/ele.14340

Agrawal A.A., Hastings A.P., Lenhart P.A., Blecher M., Duplais C., Petschenka G., Hawlena D., Wagschal V., and Dobler S., 2024, Convergence and divergence among herbivorous insects specialized on toxic plants: revealing syndromes among the cardenolide feeders across the insect tree of life, The American Naturalist, 204(3): 201-220.

https://doi.org/10.1086/731277

Agrawal A.A., Hastings A.P., Patrick E.T., Knight A.C., and Dreisbach T.A., 2025, Metabolomic variation in milkweeds shapes herbivore defense, performance, and parasite resistance, Ecology Letters.

Agrawal A.A., Salminen J.P., and Fishbein M., 2012, The coevolution of plant defenses and insect herbivores, Ecology, 93(3): 1-10.

Altizer S., and Oberhauser K., 1999, Effects of the protozoan parasite Ophryocystis elektroscirrha on the survival and reproduction of monarch butterflies, Journal of Invertebrate Pathology, 74(1): 61-68.

https://doi.org/10.1006/jipa.1999.4853

Altizer S., Hobson K.A., Davis A.K., De Roode J.C., and Wassenaar L.I., 2015, Do healthy monarchs migrate farther? Tracking natal origins of parasitized vs. uninfected monarch butterflies overwintering in Mexico, PLOS ONE, 10(11): e0141371.

https://doi.org/10.1371/journal.pone.0141371

Ben Chehida Y., van der Heijden E.S., Page E.J., Salazar C.P.A., Rosser N., Gabriela Gavilanes Córdova K., et al., 2025, Genetic parallelism underpins convergent mimicry coloration across Lepidoptera, bioRxiv, 2025-06.

https://doi.org/10.1101/2025.06.26.661542

Betz A., Höglinger B., Walker F., and Petschenka G., 2025, Regionality and temporal dynamics of sequestration and relocation of cardenolides in the monarch butterfly, Danaus plexippus, Journal of Chemical Ecology, 51(1): 19.

https://doi.org/10.1007/s10886-025-01572-8

Brower L.P., 1995, Understanding and misunderstanding the migration of the monarch butterfly (Nymphalidae) in North America: 1857-1995, Journal of the Lepidopterists' Society, 49: 304-385.

Dale J., and Stumpe M., 2014, Host plant toxicity, predator-prey interactions, and the evolution of defense: a synopsis, Ecology Letters, 17(7): 831-842.

Dalla S., and Edger P.P., 2014, Evolution of the plant Na+/K+-ATPase gene family: evidence for functional diversification, BMC Genomics, 15: 1-13.

Davey J.W., Chouteau M., Barker S.L., Maroja L., Baxter S.W., Simpson F., Merrill R.M., Joron M., Mallet J., and Dasmahapatra K.K., et al., 2016, Major improvements to the Heliconius melpomene genome assembly reveal multiple chromosomal fusions in Lepidoptera, G3: Genes, Genomes, Genetics, 6(3): 695-708.

https://doi.org/10.1534/g3.115.023655

de Roode J.C., Pedersen A.B., Hunter M.D., and Altizer S., 2008, Host plant species affects virulence in monarch butterfly parasites, Journal of Animal Ecology, 77(2): 120-127.

https://doi.org/10.1111/j.1365-2656.2007.01305.x

De-Kayne R., Gordon I.J., Terblanche R.F., Collins S., Saitoti Omufwoko K., Martins D.J., and Martin S.H., 2025, Incomplete recombination suppression fuels extensive haplotype diversity in a butterfly colour pattern supergene, PLoS Biology, 23(2): e3003043.

https://doi.org/10.1371/journal.pbio.3003043

Di Cristina G., Dirksen E., Altenhein B., Büschges A., and Korsching S.I., 2025, Pioneering genome editing in parthenogenetic stick insects: CRISPR/Cas9-mediated gene knockout in Medauroidea extradentata, Scientific Reports, 15(1): 2584.

https://doi.org/10.1038/s41598-025-85911-5

Dreisbach D., Bhandari D.R., Betz A., Tenbusch L., Vilcinskas A., Spengler B., and Petschenka G., 2023, Spatial metabolomics reveal divergent cardenolide processing in the monarch (Danaus plexippus) and the common crow butterfly (Euploea core), Molecular Ecology Resources, 23(6): 1195-1210.

https://doi.org/10.1111/1755-0998.13786

Dreisbach T.A., Hastings A.P., and Agrawal A.A., 2023, Untargeted metabolomics reveals chemical diversity and inducibility of cardenolides across milkweed species, Journal of Chemical Ecology, 49(4): 345-360.

Erickson E., Jason C., Machiorlete H., de la Espriella L., Crone E.E., and Schultz C.B., 2023, Using community science to map western monarch butterflies (Danaus plexippus) in spring, Ecology and Evolution, 13(12): e10766.

https://doi.org/10.1002/ece3.10766

Freedman M.G., and Kronforst M.R., 2023, Genomic insights into the evolution and loss of migration in monarch butterflies, Molecular Ecology, 32(8): 1950-1966.

Freedman M.G., and Kronforst M.R., 2023, Migration genetics take flight: genetic and genomic insights into monarch butterfly migration, Current Opinion in Insect Science, 59: 101079.

https://doi.org/10.1016/j.cois.2023.101079

Green D.A., and Kronforst M.R., 2019, Monarch butterflies use an environmentally sensitive, internal timer to control overwintering dynamics, Molecular Ecology, 28(16): 3642-3655.

https://doi.org/10.1111/mec.15178

Guerra P.A., Merlin C., Gegear R.J., and Reppert S.M., 2012, Discordant timing between antennae disrupts sun compass orientation in migratory monarch butterflies, Nature Communications, 3: 958.

https://doi.org/10.1038/ncomms1965

Hammer T.J., McMillan W.O., and Fierer N., 2014, Metamorphosis of a butterfly-associated bacterial community, PLOS ONE, 9(1): e86995.

https://doi.org/10.1371/journal.pone.0086995

Hemstrom W., Freedman M., Zalucki M.P., and Miller M., 2025, Novel genetic association with migratory diapause in Australian monarch butterflies, BMC Ecology and Evolution, 25(1): 43.

https://doi.org/10.1186/s12862-025-02384-w

Hoogshagen M., Hastings A.P., Chavez J., Duckett M., Pettit R., Pahnke A.P., Agrawal A.A., and de Roode J.C., 2024, Mixtures of milkweed cardenolides protect monarch butterflies against parasites, Journal of Chemical Ecology, 50(1): 52-62.

https://doi.org/10.1007/s10886-023-01461-y

Höök L., Vila R., Wiklund C., and Backström N., 2024, Temporal dynamics of faster neo-Z evolution in butterflies, Evolution, 78(9): 1554-1567.

https://doi.org/10.1093/evolut/qpae082

Li H., Handsaker B., Wysoker A., Fennell T., Ruan J., Homer N., Marth G., Abecasis G., Durbin R., and 1000 Genome Project Data Processing Subgroup, 2009, The sequence alignment/map format and SAMtools, Bioinformatics, 25(16): 2078-2079.

https://doi.org/10.1093/bioinformatics/btp352

Livraghi L., Hanly J.J., Loh L.S., Henry A., Keck C., Shirey V.M., Tsai C.C., Yu N.F., Belleghem S.M.V., and Roberts W.M., et al., 2025, Genetic basis of an adaptive polymorphism controlling butterfly silver iridescence, Current Biology, 35(9): 2154-2163.

https://doi.org/10.1016/j.cub.2025.03.028

López-Goldar X., Zhang X., Hastings A.P., Duplais C., and Agrawal A.A., 2024, Plant chemical diversity enhances defense against herbivory, Proceedings of the National Academy of Sciences, 121(51): e2417524121.

https://doi.org/10.1073/pnas.2417524121

Malcolm S.B., and Brower L.P., 1989, Evolutionary and ecological implications of cardenolide sequestration in the monarch butterfly, Experientia, 45(3): 284-295.

https://doi.org/10.1007/BF01951814

Markert M.J., Zhang L., and Reed R.D., 2016, TALEN and CRISPR/Cas9 genome editing in the monarch butterfly, Genome Research, 26(5): 792-801.

Martin A., Papa R., Nadeau N.J., Hill R.I., Counterman B.A., Halder G., Jiggins C.D., Kronforst M.R., Long A.D., and McMillan W.O., et al., 2012, Diversification of complex butterfly wing patterns by repeated regulatory evolution of a Wnt ligand, Proceedings of the National Academy of Sciences, 109(31): 12632-12637.

https://doi.org/10.1073/pnas.1204800109

Merlin C., Gegear R.J., and Reppert S.M., 2009, Antennal circadian clocks coordinate sun compass orientation in migratory monarch butterflies, Science, 325(5948): 1700-1704.

https://doi.org/10.1126/science.1176221

MonarchBase Team, 2012, MonarchBase: the monarch butterfly genome database, Database, 2012: bas028.

https://doi.org/10.1093/database/bas028

Mongue A.J., Hansen M.E.B., Gu L., Sorenson C.E., Walters J.R., and Kronforst M.R., 2017, A neo-sex chromosome in the monarch butterfly, Danaus plexippus, G3: Genes, Genomes, Genetics, 7(10): 3281-3294.

https://doi.org/10.1534/g3.117.300187

Mongue A.J., Kuperus P., and Groot A.T., 2025, Genome sequencing confirms cryptic diversity and potentially adaptive changes in gene content of Lepidoptera-infecting Apicomplexa, bioRxiv, 2025-05.

https://doi.org/10.1101/2025.05.29.656805

Mora P., Hospodářská M., Voleníková A.C., Koutecký P., Štundlová J., Dalíková M., Walters J.R., and Nguyen P., 2024, Sex-biased gene content is associated with sex chromosome turnover in Danaini butterflies, Molecular Ecology, 33(24): e17256.

https://doi.org/10.1111/mec.17256

Müller-Theissen M.L., Gottdenker N.L., and Altizer S.M., 2025, Resistance and tolerance to imperfectly specialized parasites: milkweed butterflies and their protozoan parasites, Ecology and Evolution, 15(3): e70979.

https://doi.org/10.1002/ece3.70979

Oberhauser K.S., Prysby M.D., Mattila H.R., Stanley-Horn D.E., Sears M.K., Dively G., Olson E., Pleasants J.M., Lam W.F., and Hellmich R.L., 2001, Temporal and spatial overlap between monarch larvae and corn pollen, Proceedings of the National Academy of Sciences, 98(21): 11913-11918.

https://doi.org/10.1073/pnas.211234298

Petschenka G., and Agrawal A.A., 2015, Milkweed butterfly resistance to plant toxins is linked to sequestration, not coping with toxicity, Proceedings of the Royal Society B: Biological Sciences, 282(1802): 20151865.

https://doi.org/10.1098/rspb.2015.1865

Petschenka G., Fandrich S., Sander N., Wagschal V., Boppré M., and Dobler S., 2013, Stepwise evolution of resistance to toxic cardenolides via genetic substitutions in the Na⁺/K⁺-ATPase of insects, Proceedings of the National Academy of Sciences, 110(27): 10906-10911.

Rendón-Salinas E., Alonso A., García-Serrano E., Valera-Bermejo A., and Quesada M., 2023, The monarch butterfly in Mexico: a conservation model, Current Opinion in Insect Science, 60: 101112.

https://doi.org/10.1016/j.cois.2023.101112

Sanaei E., Chavez J., Harris E.V., Alcaide T.Y., Baffour-Addo K., Bugay M.J., Adams K.L., Zelaya A., and de Roode J.C., et al., 2024, Microbiome analysis of monarch butterflies reveals effects of development and diet, FEMS Microbiology Ecology, 100(12): fiae143.

https://doi.org/10.1093/femsec/fiae143

Satterfield D.A., Maerz J.C., and Altizer S., 2015, Loss of migratory behaviour increases infection risk for a butterfly host, Proceedings of the Royal Society B: Biological Sciences, 282(1801): 20141734.

https://doi.org/10.1098/rspb.2014.1734

Semmens B.X., Semmens D.J., Thogmartin W.E., Wiederholt R., López-Hoffman L., Diffendorfer J.E., and Pleasants J.M., 2016, Quasi-extinction risk and population targets for the eastern migratory monarch butterfly (Danaus plexippus), Scientific Reports, 6: 23265.

https://doi.org/10.1038/srep23265

Shen K., 2024, The use of CRISPR-Cas9 reveals the major genes that control wing patterning and coloration in Nymphalidae butterflies.

Thogmartin W.E., Wiederholt R., Oberhauser K., Drum R.G., Diffendorfer J.E., Altizer S., Taylor O.R., Pleasants J., Semmens D., and Semmens B., et al., 2017, Monarch butterfly population decline in North America: identifying the threatening processes, Royal Society Open Science, 4(9): 170760.

https://doi.org/10.1098/rsos.170760

U.S. Fish and Wildlife Service, 2020, Species status assessment report for the monarch butterfly (Danaus plexippus), U.S. Department of the Interior.

Uffelmann E., Huang Q.Q., Munung N.S., De Vries J., Okada Y., Martin A.R., Martin H.C., Lappalainen T., and Posthuma D., 2021, Genome-wide association studies, Nature Reviews Methods Primers, 1(1): 59.

https://doi.org/10.1038/s43586-021-00056-9

Zhan S., Merlin C., Boore J.L., and Reppert S.M., 2011, The monarch butterfly genome yields insights into long-distance migration, Cell, 147(5): 1171-1185.

https://doi.org/10.1016/j.cell.2011.09.052

Zhan S., and Reppert S.M., 2013, MonarchBase: a genomic resource for the monarch butterfly, Nucleic Acids Research, 41(D1):

https://doi.org/10.1093/nar/gks1057

Zhan S., Reppert S.M., and Kronforst M.R., 2020, Chromosome-level genome assembly of the monarch butterfly (Danaus plexippus), GigaScience, 9(5): giaa042.

Zhan S., Zhang W., Niitepold K., Hsu J., Haeger J.F., Zalucki M.P., Altizer S., de Roode J.C., Reppert S.M., and Kronforst M.R., 2014, The genetics of monarch butterfly migration and warning colouration, Nature, 514(7522): 317-321.

https://doi.org/10.1038/nature13812

Zhang L., and Reed R.D., 2016, Genome editing in butterflies reveals that spalt promotes and Distal-less represses eyespot colour patterns, Nature Communications, 7: 11769.

https://doi.org/10.1038/ncomms11769

Box 1 Conservation Implications of Monarch Genomics and Metabolomics

Recent advances in monarch genomics and metabolomics provide an opportunity to move beyond abundance-based conservation toward strategies that preserve functional, adaptive capacity across the migratory range.

1. Cheically informed habitat restoration. Milkweed species and populations differ substantially in cardenolide composition, inducibility, and secondary metabolite diversity, with direct consequences for monarch survival, chemical defense, and resistance to the protozoan parasite Ophryocystis elektroscirrha (OE) (de Roode et al., 2008; Agrawal et al., 2012; Petschenka and Agrawal, 2015). Metabolomic screening of milkweeds used in restoration projects could identify plant chemotypes that optimize larval performance while enhancing parasite resistance and predator deterrence, avoiding one-size-fits-all planting strategies.

2. Preserving adaptive genetic variation for migration. Population genomic studies indicate that migration, diapause, lipid storage, and navigation are polygenic traits shaped by subtle allele-frequency shifts across many loci (Zhan et al., 2014; Freedman and Kronforst, 2023). Conservation actions that maintain connectivity among breeding, migratory, and overwintering regions are therefore essential to preserve adaptive alleles associated with circadian timing, endocrine regulation, and metabolic endurance.

3. Metaolites as early-warning indicators. Metabolomic profiles-such as cardenolide sequestration patterns, lipid reserves, and stress-response metabolites-may provide sensitive indicators of physiological condition and migratory readiness before population declines are detectable through census data alone (Semmens et al., 2016; Thogmartin et al., 2017). Integrating metabolomic monitoring into long-term surveys could improve detection of sublethal stress caused by climate extremes, pesticide exposure, or host-plant mismatch.

4. Integrating host–parasite–microbiome dynamics. Chemical defense, parasite resistance, and gut microbiome composition interact to shape monarch fitness. Restoration strategies that consider host-plant chemistry alongside disease pressure and microbial interactions may reduce parasite prevalence and improve survival during migration and overwintering (de Roode et al., 2008; Hammer et al., 2014).

5. From molecular insight to management practice. By integrating causal genomics, pan-genomics, regulatory neurogenomics, and metabolomics, conservation efforts can prioritize not only habitat quantity but also genetic, chemical, and physiological quality. Such mechanistically informed strategies are likely to be more robust to environmental change, supporting long-term persistence of migratory monarch populations in an increasingly variable landscape (U.S. Fish and Wildlife Service, 2020; Erickson et al., 2023).